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Abstract 

The non-factorizable terms in the operator product expansion have been recognized 
as one of theoretical obstacles for precision determination of the mixing parameter V u b 
from semileptonic B decays. It is pointed out here that the recent CLEO data on the 
parameters of the heavy quark expansion Ai and A, combined with a theoretical bound 
on Ai strongly favor the existence of a sizeable contribution of non-factorizable terms 
in semileptonic decays of D mesons. Thus these terms are likely to solve the long- 
standing problem of the deficit of semileptonic decay rate of the D mesons, and with 
better data their magnitude can be determined and used in studies of the parameter 

Vub- 



The determination of the weak mixing parameter V u b from the data on inclusive semilep- 
tonic decays of B mesons B — > X u t v requires kinematical cuts in order to discriminate 
against the dominant contribution of B — > X c lv (see e.g. in the review talks ]]], 0)- The 
inclusive decay rate in the kinematically restricted part of the phase space however is sen- 
sitive to effects that being relatively small in the total rate are dominantly concentrated in 
the kinematical region of interest. A typical example of such behavior is presented by the 
method P, f§] using a cut in the invariant mass of the lepton pair: q 2 > (Mb — Mp) 2 , which 
leaves only about 20% of the total inclusive decay rate of B — > X u lv. When considered 
within the operator product expansion this inclusive decay rate however receives a contri- 
bution from a matrix element over the B mesons of a non-factorizable four-quark operator. 
More specifically, this contribution reads as (see e.g. || || and references therein) 

iZ 7T 

where fs is the B meson annihilation constant, and the phenomenological parameters "bag 
constants" Bi and B 2 parameterize the matrix elements of four quark operators over the B 
mesonQ 

(51(^(1-75)^(^(1-75)6)15) = &p-B t , 

(B\b(l- l5 )u)(u(l+ l5 )b)\B) = f -^B 2 . (2) 

In the factorization limit, where the matrix elements are saturated by the vacuum insertion 
between the two-quark factors, both constants B are equal to one and cancel in the expression 
([]]). However the deviations from factorization are expected to be at a ~ 10% level. (E.g. 
these deviations are generally 0(N^ 2 ) in the limit of large number of colors Nc-) 

At the expected level of about 10% of breaking of the factorization relations, the natural 
magnitude of this contribution amounts to only about 2.5% of the total rate, but it is all 
concentrated near the maximal value of q 2 , formally 5(q 2 — m 2 ), and thus is clearly greatly 
enhanced relative to the rate in the restricted part of the phase space. For this reason a 
better understanding of the magnitude of the non-factorizable terms is of a great importance 
for an accurate determination of \V u b\ from the data. 

The differences of the non-factorizable terms between heavy mesons with different light 

quark flavors, i.e. the flavor non-singlet part of these terms, can be tested |J by measuring 

lr The nonrelativistic normalization: e.g. (B\b'b\B) = 1, for heavy quark operators is used throughout this 
paper. 



1 



the differences of the semileptonic decay rates either between the B mesons or between 
the D mesons. An evaluation || of the non-singlet part of a similar non-factorizable term 
(although with a somewhat different color structure) from the difference of lifetimes of D s 
and D mesons confirms that the deviations from factorization are indeed at a ~ 0.1 level. 

An estimate of the light-flavor singlet part however requires an evaluation of the absolute 
decay rates of heavy mesons, rather than their differences, which also contain other, formally 
leading in the heavy quark expansion terms that are flavor- independent. Naturally, extract- 
ing the contribution of the non-factorizable terms from the absolute decay rates requires a 
quite precise knowledge of the parameters involved in the leading, perturbative term, as well 
as in the 0{m^) terms. This essentially leaves the semileptonic decays of D mesons as the 
only appropriate testing ground for the flavor singlet part of the non-factorizable terms. 

The problem of explaining the total semileptonic decay rates of D mesons has been 
discussed in the literature in the past |7|, ||, where it has been noticed that the theoretical 
decay rate, scaling with mass of the charmed quark essentially as m^, is significantly lower 
than the experimental value if m c is taken to be around 1.4 GeV, as suggested by the QCD 
sum rules for charmoniumf^] and possibly by other considerations. In order to remedy 
this contradiction Chernyak |7J suggested that the appropriate value of m c is substantially 
larger, m c « 1.65 GeV, while Blok, Dikeman and Shifman ||, keeping m c not larger than 
about 1.4 GeV, considered and dismissed a possible contribution of non-factorizable terms, 
but rather attributed the contradiction to possible violation of the quark-hadron duality. 
Certainly a violation of the duality and a general failure of the OPE at a scale as low 
as the charmed quark mass is a logical possibility, in which case a consideration of the 
B — > X u i v decay rate with the cut q 2 > (Mb — Mp) 2 would also be of little meaning, since 
the scale involved in this case fi c ~ m c (1 — m c /2mt ) ) |4| is yet somewhat lower. However, 
since there are no other indications of such failure, it looks more constructive to still use 
the heavy quark expansion and to explore the possibility that it is the contribution of the 
non-factorizable terms that reconciles the experimental semileptonic widths of the D mesons 
with the theoretical expression using a moderate value of m c . The latter possibility, also 
discussed in Ref. ||, can now be addressed with better certainty, due to recent CLEO data 
[|To|1 on the parameters Ai and A of the heavy quark expansion. 

One can notice that when applied to semileptonic decay rates of the D mesons, similarly 
to eq.(|l|) the contribution of the non-factorizable terms can be written as 

8V sl {D) = ^ ; f o m>D (\V CS \ 2 5B S + \V cd \ 2 5B d ) , (3) 
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where 5B q is an analog of the difference B2 — B\ for operators involving a light quark q, 
instead of the u quark in eq.(fj). In the leading order in heavy quark mass the parameter 
SB q can be written in terms of the matrix element of an operator involving only spatial 
components (in the rest frame of the heavy meson) of the V — A currents: 



2{P Q \(Q^(l- l5 )q)(m- l5 )Q)\P Q ) = f 2 P m P 6B q , (4) 

with Pq standing for a pseudoscalar meson containing a heavy quark Q. Numerically, the 
expression (131) gives 



m c \ 2 ( f D \ 2 (\V cs \ 2 5B s + \V cd \ 2 5B, 



sr,(o) = <>.«.»»-' [^) p' ""' ) ■ (5) 

Thus with 'natural' values 5B s ,5Ba ~ 0.1 the effect of the non-factorizable terms easily 
reaches about one half of the experimental semileptonic decay rate, e.g. T s i(D°) = 0.164 ± 
0.007 ps^ 1 ||1 1|| . Therefore an analysis of these rates necessarily should include the non- 



factorizable terms even at their expected suppressed level. 

As already mentioned the estimate of about 1.4 GeV for the 'pole' mass of the charmed 
quark originates from the charmonium sum rules || . Subsequent development of this method 
had however revealed that the perturbative expansion in a s for the propagation of the cc 
pair, used in derivation of the sum rules, involves unusually large coefficients in the order a 2 s 



14 . This might cast a general doubt on reliability of extraction of the parameter m c from 
the sum rules. 

An alternative way of extracting the heavy quark mass parameter from the data is based 
on the mass formula for a heavy pseudoscalar meson: 

M P = m Q + A + ^i + 0(m- Q 2 ) , (6) 
2m Q v 

where A is the mass shift due to light degrees of freedom in the meson, and fi 2 and fi 2 
describe the kinetic and chromomagnetic energy of the heavy quark inside the meson: fi 2 = 
-(P\(QD 2 Q)\P), fj% = (P\{Q(a ■ B)Q)\P) with L> M being the QCD covariant derivative and 
B being the chromagnetic field operator (see e.g. in the review |13[)- The parameter fi 2 is 
readily found from the mass splitting of vector and pseudoscalar mesons, n 2 ~ 0.37GeV 2 , 
while the values of the parameters A and [i 2 are less certain. 
The value of fi 2 is however constrained by the inequality ]14 



&>fi, (7) 
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which follows, in terms of the nonrelativistic expansion for a heavy quark, from the non- 
negativity of the Pauli Hamiltonian —{a ■ D) 2 = —D 2 — (a • B), or from a more general 
inequality, not necessarily related to the heavy quark expansion, 

= (P\(QQ)\P) < (p\{Q ] Q)\p) = i • (8) 

When applied to the expansion (El) the latter general inequality reduces to the constraint 



Recently CLEO has published |T(| an analysis of the parameters A and Ai = —fi 2 from 



their data on the inclusive photon spectrum in the decays B — > X s 7 and on the second 
moment of the hadronic invariant mass in the semileptonic decays of B mesons^. The results 
of this analysis can be converted, using the mass formula (|6|), into a value of m c , which in 
turn can be used in the expression for the total semileptonic decay rate of the D mesons and 
thus enables an evaluation of the contribution of the non-factorizable terms from comparison 
with the experimental rate. A calculation along these lines is to be presented in the rest 
of this paper. The additional element in this analysis as compared to Ref.Jl^] is that the 
constraint (^), Ai < —0.37, is explicitly taken into account. 

The expression for the semileptonic decay rate of a D meson, including the perturbative 
QCD corrections up to two loops (16| and the terms in heavy quark expansion of order m~ 2 
fl7|l as well as the non-factorizable terms of order m~ 3 can be written as 

| V C S\ I 1 - ° n I T I V cd \ A 

V m cJ 

l-^]+5T sl (D), (9) 

where a s = a s (m c ), 6T s[ (D) is given by eq.(|3]), and a certain inaccuracy has to be admitted 

in the treatment of the cross terms between e.g. the radiative corrections and the effect of 

the finite mass m s of the strange quark or between the radiative corrections and a part of 

the 0(m~ 2 ) terms. This inaccuracy however is at the level of other uncertainties involved in 

eq.(P), e.g. due to higher perturbative terms, or the experimental uncertainties in the data, 

and can be safely neglected in the present calculation. 

2 It can be noticed that the general inequality (]|) implies the bound |l5| on the difference of the b and c 
quark masses: mj, — m c > Mb — Md ~ 3.41 GeV, which does not rely on the heavy quark expansion. 

3 It should be noted, for the sake of rigor, that depending on definition of perturbative renormalization 

contributions there may arise a difference between Ai and — fi^, as discussed in Ref.|13|. However in the 




context of the analysis of Rcf.(T^], no such difference arises, and Ai and — /x 2 are equivalent. We switch to 
the notation Ai in the rest of this paper in order to facilitate following the relation to Ref . . 
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Figure 1: The lines of constant relative contribution 5 of non-factorizable terms to T s i(D) 
overlayed on the CLEO results |L0| in the physical region Ai < —0.37 at a s {m c ) = 0.25 and 
cn s (r?T c ) = 0.3. The elliptical arcs correspond to the indicated levels of x 2 m the CLEO data. 
The filled circle is centered at the minimum of x 2 of about 1.2 in the physical region. The 
lines in the left part of the plots correspond to the indicated values of 5. 



The equation (H) can be used for evaluating the contribution of the non-factorizable terms 
as a function of Ai and A, since these also determine the mass m c through the mass formula 
and the actual rate T s i(D) is known experimentally. The results of such evaluation 
are shown together with the CLEO data |T(J in Fig.l in terms of the dimensionless ratio 
5 = 5T s i(D)/T s i(D ), where the experimental value is used for T s i(D°) in the denominator. 
Only the physical region Ai < —0.37 is shown in the plots. Although the central point of the 
CLEO result lies above this region, the minimal value of x 2 m the physical region is about 1.2, 
so that the data are well consistent with the constraint (0). The parameters corresponding to 
this point of maximum likelihood are A = 0.456 GeV and (naturally) Ai = —0.37, saturating 
the bound (^). According to the mass formula (§) these values translate into m c m 1.41 GeV 
and, when used in eq.(^) for the T s i(D), the maximum probability corresponds to S ~ 0.5 if 
<^ s { m c) — 0.25 and to S ~ 0.6 if a s (m c ) = 0.3. (Clearly, the estimated value of S grows with 
a s (m c ), since the QCD radiative corrections suppress the perturbative and 0{m~ 2 ) terms 
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in the rate.) It can be also noted that the point of maximum likelihood corresponds to the 
'best value' of the appropriate 'pole' mass of the b quark pa 4.82 GeV, in a reasonable 



agreement with the evaluation from the QCD sum rules for T resonances fT8|j . 

The most likely values of m c and 5 suggested by the CLEO data are in a very good 
agreement with theoretical expectations from other considerations. In particular the value 
of S around 0.5 — 0.6 agrees with the expectation that the suppression of the non-factorizable 
terms is about 0.1. As seen from the plots in Fig.l the value 5 = 0, and hence the mass m c 
around 1.6 — 1.65 GeV, has quite low likelihood and is only marginally consistent, at a level of 
about 3a, with the CLEO data and the constraint (0). Clearly, a more statistically significant 
evaluation of the parameters of the heavy quark expansion and thus of the quantity 5 requires 
improvement of both the theoretical and experimental accuracy in an analysis along the lines 
of that in Ref.[|T0f supplemented by the remarks presented in this paper. In any event, a 
precision determination of the weak mixing parameter \V u b\ from a constrained spectrum of 
semileptonic B decays is quite unlikely before the issues related to semileptonic decay rates 
of D mesons are resolved both theoretically and experimentally. 

I thank A. Czarnecki for an update on the status of perturbative calculations of semilep- 
tonic decay rates, R. Poling for a discussion of the CLEO work |TU[ and D. Cronin-Hennessy 
for e-mail communication concerning the results of the CLEO analysis and the parameters 
of the final plot in their paper. This work is supported in part by DOE under the grant 
number DE-FG02-94ER40823. 
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